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Abstract. The extinction laws at 7/xm and 15/im are derived for more than 120 
sightlines in the inner Galactic plane based on the ISOGAL survey data and the 
near-infrared data from DENIS and 2MASS. The tracers are the ISOGAL point 
sources with [7]-[15]<0.4 which are RGB tip stars or early AGB stars with moder- 
ate mass loss. They have well-defined intrinsic color indices (J-Ks)o, (Ks-[7])o and 
(Ks-[15])o- By a linear fitting of the observed color indices Ks-[7] and Ks-[15] to 
the observed J-Ks, we obtain the ratio between the E(Ks-[7]) and E(Ks-[15]) color 
S-H ' excesses and E(J-Ks). We infer the selective extinctions at 7 and 15pim in terms of 



the near-infrared extinction in the Ks band. The distribution of the derived extinc- 
tions around 7 micron (A7) is well represented by a Gaussian function, with the 
peak at about 0.47Ak s and ranging from 0.33 to 0.55Ak s (using the near-infrared 
extinctions of Rieke & Lebovsky 1985). There is some evidence that A7/AKS may 
vary significantly depending on the line of sight. The derived selective extinction 
at 15/xm suffers uncertainty mainly from the dispersion in the intrinsic color in- 
dex (Ks-[15])o which is affected by dust emission from mass-losing AGB stars. 
The peak value of A15 is around 0.40Ak s . 

Key words, stars: AGB and post-AGB -stars: late-type - stars: mass loss - ISM: 
extinction - Infrared: stars 



1. Introduction 



In the mid-infrared waveleng th range, th e interstellar extinction is basically composed 
of two components [see e.g. iDraind l|1989|) ] . One is the continuum extinction decreasing 
Send offprint requests to: Biwei Jiang (bjiang@bnu.edu.cn) 
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with wavelength which can be attributed to graphite or PAH carbonaceous dust grains; 
the other relates to the spectral features around 10/im and 18/mi due to silicate dust 
grains. The extinction around 7/zm is rather controversial, partly due to the fact that 
this wavelength is where the decreasing continuum extinction meets the rising silicate 
extinction peak. Before ISO (Infrared Space Observatory) was launched, there were sev- 
eral measurements of near-infrared and mid-infrared extinctions which could be used to 



derive that at 7/xm. 



Rieke fc Lebofskvl 1|1985|) observed o Sco and a nu mber of stars in 



the Galactic cent er and derived t he ex tinction law in the 1-13/zm range; 



reobserved o Sco: 



Landi ni et al 



Whittejj JL 



1984) observ ed hydrogen recombination lines from the 



obscured HII region G333.6-0.2.|DraineJ l)l989l) concluded that the available observational 
evidence was consistent with the power law A\ oc A~ 175 out to A w 7/xm (which means 
the extinction at 7jum is about 0.13Ak) and that the various extinction curves from 
different observations appear to be essentially universal for wavelengths A > 0.7/xm. 

After the launch of ISO, new determinations were made of the extinction around 
7/xm. Using ISO observations of hydrogen recombination lines, 



Lntzl U. 




Lutz et al 



and 



found that the extinction toward Sgr A* does not decline with increasing 
A in the 4-8/^m region, and obtained an extinction of about 0.45Ak at 6.8/xm. However, 
from their studies u sing H2 rovibrational lines in the Orion molecular cloud outflow, 



Rosenthal et al 



found that the extinction continues to decline with increasing 
A to a minimum at about 6.5/^m, and the extinction at 7/Ltm is about 0.15Ak, though 
their lower value may be attributed to the relatively modest extinction of the field they 
analyzed, Ak s being only 1 mag. Moreover, by comparing the intensity of infrared dark 
clouds and their neighboring emission and assuming a simple radiative transfer model, 
( 200 11) derived an average o pacit y radio of 7/xm to 15Mm o f 0.7±0.1. The 



Hennebelle et al 



later two results from 
with the re sult by 
results from 



Rosenthal 



et al 



( 2000) and 



Hennebelle et al 



12001) are consistent 



Dramel 1 1989 ). The issues of whether th e discrepancy between t he 



Lutz (1999) and from 



Rosenthal et al 



(200QJ) or 



Hennebelle et al 



(20011) is 



real or not and whether it could result from dust properties that depend on the line of 
sight need more investigation. 

With data a vailable from numerous lines of sight in the Galactic disk, the ISOGAL 



and GLIMPSE l)Beniamin et all 2003: 



Churchwell et al . 2004) projects provide an oppor- 



tunity to study the dependence of the extinction on the interstellar dust environments. 

The ISOGAL project is a survey of 263 fields, covering about 16 square degree area 
in total, in the Galactic pla ne. It made use of th e ISOCAM camera aboard the ISO space 
telescope at 7/im and 15/im l Omont et al.12003 ). Furthermore, most of the ISOGAL fields 



in th e southern hemisphere were also observed by the DENIS project l)Epchtein et al 



1999J) . in the IJKs bands. The ISOGAL data w ere systematically merged with DENIS in 



a five-wavelength catalog 



In the meantime, another ground-based 



near-infrared sky survey project, 2MASS, published data covering almost the whole sky, 
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including the ISOGAL sky area. In combination with the DENIS and the 2MASS data, 
ISOGAL has great potential for studying the extinctions at 7 and 15/im in the Galactic 
plane. 

71 colors obtained from D ENIS and ISOGAL 



By combining J 


Omont et al. 


2003 





Schullcr et al 



2003), we IjJiang et al 



2003) derived the selective 



extinction at 7 and 15/im along the line of sight towards an ISOGAL field at /=18.3 and 
6=-0.35, obtaining A7 ~ 0.4Ak s - Althou gh this res ult suffers from some u ncertainty, it 



is more consistent with 



Lutz et al 



(1996) than with 



Rosenthal et al 



(2000). 



The most recent result concerning the mid-infrared extinction is from Spitzer. The 
GLIMPSE program observed a lar g e part of the Galac t ic disk in four bands between 3.6 

1 2005) found, along two lines of 



and 8/im Ijlndebetouw et al 



2005). 



Indebetouw et al 



sight , that the extinction across the 3-8/im wavele ngth range is flat, which agrees with 



Lutz (1 



as well as ourselves 



Jiang et al 



2003). 



Given that the wide areas covered by ISOGAL and, more systematically, by 
GLIMPSE, include various sightlines (and thus various extinction environments) one 
can check whether the mid-IR extinction varies. The present paper aims at checking the 



Jiang et al 



extinction at 7 and 15 /mi. We have used in principle the same method as in 
I 2003) to determine the extinction law at 7 and 15/im from the ISOGAL fields. Since the 
ISOGAL fields are within the Galactic plane and lie mostly in the inner disk, this work 
mainly relates to the central parts of the Galactic plane. 



2. Database 



2.1. The ISOGAL PSC data 

This work is base d on the ISOGAL-DENIS Point Source Catalog Version 1.0 



( Schuller et al 



2003) which includes more than 100,000 point sources detected at 7/im 
and/or 15/im in the ISOGAL survey. In addition to the standard data reduction by the 



CIA package (IQtt et al 



1997J) of the ISOCAM images, the point sources are extracted 
by a dedicated PSF fitting procedure. Because most of the ISOGAL survey areas are 
near the Galactic plane < 1°) and have high source density, the source extraction is 
confusion limited and disturbed by strong background radiation. The photometric error 
is usually better than 0.2mag and rises above this value at the ISOGAL sensitivity limit 
of about lOmJy. The ISOGAL PSC catalog consists of three parts. The first part includes 
163 fields observed at both 7/im and 15/im, with fields denoted by "FC". The second 
part includes 43 fields observed only at 7/im, denoted by "FA" . The third part includes 
57 fields observed only at 15/im, denoted by "FB" . 

Unfortunately, the ISOGAL observations did not use the same filters through all 
the fields: towards about 75% of them the broad-band filters LW2 and LW3 were used, 
while towards the other 25%, the LW5, LW6 and LW9 narrow-band filters were used 
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(Schullcr et al 



2003). The reference wavelengths of the LW2 and LW3 filters are 6.7/im 
with a bandwidth of 3.5/im and 14.3/mi with a bandwidth of 6.0/im, respectively. The 
LW5, LW6 and LW9 filters have reference wavelengths at 6.8/zm, 7.7/zm and 14.9/i m, with 



2003). Thus, 



bandwidths of 0.5/im, 1.5/xm and 2.0/im, respectively ijBlommaert et al. 
what we call "7/im" or "15/im" is an approximation and should not be confused with 
the precise central wavelengths of the bands. Detailed qualitative information on the 
ISOGAL filters, including the reference wavelengths, bandwidths, zero poi nt magnitudes 



and flux densities, and total observed area, can be found in Table 1 of 



(2003). 



Schuller et al 



2.2. The ISOGAL-DENIS PSC data 

The nominal position error is about 0.5" in the DENIS survey, while it is several arc 
seconds in the ISOGAL survey. Cross-identification with DENIS has thus improved signif- 
icantly the positional accuracy of the ISOGAL point sources. The final position accuracy 
of the ISOGAL-DENIS PSC is normally about 1 arc second. To deal with the confusion 
problem, the Ks limiting magnitude was cut to give a source density of 72000/deg 2 for the 
ISOGAL 3"pfov [pixel field-of-view] observations or 18000/deg 2 for the IS OGAL 6"pfov 



observ ations. The details of the ISOGAL-DENIS catalog are described in 



Schuller et al 



(2003). Although there are over 100,000 sources in the ISOGAL-DENIS PSC, only data 
of high quality are used in the present work. We chose only sources whose association 
flags are 4 or 5, implying that the cross identifications are reliable. When there are 
two sources within the search radius, the closest DENIS source is associated with the 
ISOGAL source, and the association flag is usually reduced by 1 (from 5 as the high- 
est quality). Such cases are mostly dropped from the following analysis. In addition, a 
constraint was also applied to the photometric quality of the DENIS PSC sources. The 
correlation factors in both J and K bands were required to be larger than 0.85. Since 
DENIS surveyed only the southern sky, it did not cover all the ISOGAL fields. The final 
number of retained associated sources is 50724 distributed in 220 fields, among which are 
8552 sources in 42 FA fields (7/um only), 7107 sources in 56 FB fields (15/im only) and 
35065 sources in 122 FC fields (both 7 and 15/xm). 

2.3. The ISOGAL-2MASS PSC data 



The Two-Micron All Sky Survey, ! 2MASS. lSkrutskie et all |l997L is similar to DENIS in 



that it is also a ground-based near-infrared survey. However, 2MASS differs from DENIS, 
mainly in that it covers the whole sky including all the ISOGAL fields. In addition, 
2MASS observed in the JHKs bands, i.e. including the H band but not the I as in 
DENIS. The depth of the 2MASS photometry is slightly greater than that of the DENIS 
photometry; nominally the limiting magnitude for point sources is Ks=14.3mag compared 
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to 14 mag in DENIS. As the field of view per pixel is similar, the positional accuracy of 
the 2MASS PSC is of the same order as that of D ENIS. 

The 2MASS PSC database was released in 2003 ( Cutri et al.l2003 ) , after the ISOGAL- 



DENIS PSC. Th e cross-identification p rocedure was the same as for the ISOGAL-DENIS 



case, described in 



Sch uller et al 



I 2003). There are three major steps. First, we determined 
the offsets of every ISOGAL field according to the positions provided by the 2MASS 
PSC. As with all the southern ISOGAL fields, this correction was accomplished by cross- 
identification with the DENIS PSC. As the 2MASS astrometric accuracy is comparable 
to that of DENIS (both have a nominal error of 0.5"), this would not have been any sig- 
nificant improvement to the positional accuracy. Consequently, this step is applied only 
to the 42 northern fields which were beyond the DENIS coverage. The offsets of these 
ISOGAL fields to the 2MASS PSC sources range as high as 4 to 5 arcseconds though 
usually they are less than 3", similar to the values of the ISOGAL field offsets to the 
DENIS PSC sources, as expected. After removing the offsets, the positional accuracy 
of sources in the northern ISOGAL fields is about 1.5". The second step is to search 
for the 2MASS counterpart of each ISOGAL point source. The search radius is 1.8" or 
3.5" for sources in the ISOGAL fields with pfov of 3" and 6" respectively. If there is 
more than one 2MASS source within the search radius, no 2MASS source is assigned 
to the ISOGAL source (although the closest one was taken to be the ISOGAL-DENIS 
cross-identification). Also rejected are the sources whose positional differences between 
ISOGAL and 2MASS are bigger than 3a, where a means the rms ISOGAL-2MASS posi- 
tion difference in the corresponding field (after field offset correction). The third step is 
to try to avoid the confusion problem which was met with in the cross-identification with 
the DENIS sources. It is likely to be more serious because of the slightly higher sensitivity 
of 2MASS. The 2MASS source density is limited to 72000/deg 2 or 18000/deg 2 for the 
ISOGAL fields with pfov 3"or 6"respectively, which corresponds to a formal random as- 
sociation probabi lity of 5%, implying that the number of spurious associations should be 



smaller than 1% (Schuller et al 



2003). This step is accomplished by setting the limiting 



magnitudes in the Ks band above the 2MASS limit (which varies from field to field and 
ranges from about 12 to 14 mag). The constraints on the association radii and confusion 
limits are similar to good association quality flags (4 or 5) in the ISOGAL-DENIS PSC. 
In addition to these steps, which are similar to those involved in the ISOGAL-DENIS 
association process, the 2MASS PSC sources are required to have reliable photometric 
results by only accepting photometry quality flag "AAA" in the JHKs bands, which 
means S/N > 10 in all bands. The read flags in all three bands were also required to 
be good (i.e. "Rflag" being 1, 2 or 3), which indicates best quality detection, photom- 



etry and astrometry IjCutri et al 



2003). Although these two criteria may be redundant 



in some sense, they are consistent and guarantee the quality of data. Under these con- 
straints, the resultant ISOGAL-2MASS PSC contains 50653 sources distributed in all 
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the 263 fields. These sources are composed of 8112 in 43 FA fields, 6623 in 57 FB fields 
and 35918 sources in 163 FC fields; or 39017 7/j,m sources in 205 fields and 23359 15/j,m 
sources in 218 fields. Compared to the ISOGAL-DENIS PSC sources, the number of fields 
is increased by 43 thanks to the 2MASS coverage of the northern ISOGAL fields. On 
the other hand, the total numbers of associated sources are comparable, which means 
that there are fewer sources on average in each field. This decrease in average density of 
sources in the ISOGAL-2MASS catalog can be attributed to the stricter constraint on 
photometric quality for the 2MASS PSC sources. This ISOGAL-2MASS PSC catalog is 
available via the VizieR Service at the Centre de Donnees Astronomiques de Strasbourg 
(CDS, http://vizier.u-strasbg.fr/viz-bin/VizieR/). 



3. Method 

3.1. Assumptions 

We use the method described in detail by 



Jiang et al 



I 2003J) to derive the extinction 



values. In the following sections, we adopt the following symbol convention: magnitudes 
in the J,Ks, bands and at 7 and 15/mi are denoted by J, Ks, [7] and [15] respectively; 
intrinsic color indexes m\ — m\ by C AiA2 ; observed color indexes by Ca^; the color 
excess of the observed to the intrinsic color C" iAs , — C^Aa by E(Ai — A2); the extinction 
by Aa- The method is based on two assumptions. The first assumption is that the group 
comprising RGB-tip and early AGB stars have approximately constant intrinsic color 
indexes Cj Ks , C^ s7 , C^ sl5 . The second is that the ratio of color excess E(Ks-[7]) to E(J- 
Ks),i.e. k7 =E(Ks-[7])/E(J-Ks), is constant in one field. Then the extinction at 7/im can 
be calculated by a linear fitting of the observed color indexes Ck S 7 to Cjk s , i-e. 

Cks7 - Cks7 = k 7 X (CjKs - Cj Kb ) (1) 

^Ks -A 7 = k 7 x (A] - A Ks ) (2) 
The same formula can be applied to the 15/im case. 



3.2. Selection of sources 



Two types of sources are to be excluded because their intrinsic colors differ from our 
standard red giant sources: i) dusty young stellar objects and AGB stars with thick 
circumstellar envelopes, and ii) foreground early-type stars. The former are distinguished 
by their color index Cj 15 which is usually larger than 0.4. Thanks to the small extinction 
at both 7 and 15/im, the observed co l or ind ex C715 may be used as a good approximation 



Felli et al 



1 2002) calculated the influence of the circumstellar 



to the intrinsic one Cy 15 . 
matter on C7 15 and found that the sources having this color index greater than 0.4 may 
have thick circumstellar matter and therefore redder intrinsic C^ sl5 and C^ s7 . However, 
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Fig. 1. Distribution of the ISOGAL fields in the Galactic plane. The filled circles denote 
the fields for which reliable results on the 7/xm extinction were obtained and the open 
circles denote the fields for which no reliable results on the 7/zm extinction could be 
obtained. 



selection by color is applicable only to the sources in the "FC" fields which were observed 
at both 7/im and f 5/im. It is not suited to the sources in the "FA" or "FB" fields observed 
at only one wavelength since no color index C715 can be derived. Most foreground early- 
type stars can be excluded simply by dropping the sources with Cjk s <1.2. 

Some fields lie in directions with small extinction, as shown by th eir small E(J — Kg ) 

Jiang et alJ feooa) . 



values. From the analysis of the extinction law at 7/im and 15/im by 
A7 is about 0.18E(J — Ks) and A15 is about 0.16E(J — Ks). Thus, an extinction of the 
order of E(J — Ks)= 1 mag results in A7 ~ 0.18 mag, comparable to the photometric 
error. In order to be able to obtain significant extinction in these mid-infrared bands, 
only the fields with much greater near-infrared extinctions are targeted. Only sources 
with large values of the color excess E(J — Kg) are chosen so that their extinction at 
mid- infrared wavelengths will be significant. This will be further clarified in the next 
section. 

In summary, a source to be selected should have: good photometry quality both 
in the ISOGAL mid-infrared and DENIS or 2MASS near-infrared bands, good near- 
infrared-mid-infrared association flags, and a moderate mass loss rate with C715 <0.4. In 
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addition, for a field to be qualified, there must be enough sources which suffer significant 
and detectable extinction in the mid-infrared bands. 



FC-0S992+00012 FC+0 1 694+0008 1 




-1 0123403468 10 

Ks-[7] (mag) Ks-[7] (mag) 

Fig. 2. Two example fields showing all the ISOGAL-2MASS associated sources. The 
crosses denote the sources kept for final fitting. The other symbols denote various kinds 
of sources rejected before fitting, specifically, the asterisks for J-Kg <f.2, open circles for 
[7]-[15]>0.4, and triangles for deviations from the fitted line greater than 3a. The figure 
on the left shows a good-quality case, FC-02992+000f 2, with small a^ 7 ; the figure on the 
right shows a mediocre case, FC+Of 694+00081, with relatively large error in k?. In the 
right figure, the points with deviation between 2a and 3a are shown also. The results 
from fitting by the 2a and 3a rejections (denoted by dash and solid lines respectively) 
are very close to each other fsee l3.3[) . 



3.3. Fitting Procedure 

Fitting is done in two steps. First, we make a preliminary robust fitting to all the points 
which satisfied the previous criteria. In comparison with normal least-square linear fitting, 
the robust fitting method, which minimizes the absolute deviation, is less sensitive to 
small number of large outliers and more appropriate for our cases which usually have 
fewer sources with large color indexes. This fitting distinguishes the sources with large 
deviation, i.e. deviation larger than 3a from the fitted line. These sources were rejected 
because they have little probability of obeying the linear relation from the statistical point 
of view; on the other hand, most of them have relatively large Cks7, which indicates the 
likelihood of extinction arising from circumstellar dust envelopes. The linear fitting was 
carried out both ways, one fitting Ck S 7 (as y) to Cjk s (as x) and the other fitting Cjk s 
(as y) to Cks7 (as x); a source deviating from either fitting larger than 3a was rejected. 
We also tested 2a as the critical deviation for rejection and found the results were little 
changed. As an example, the right panel in Fig. [5] shows the sources with deviation 
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between 2a and 3cr as well as the results by fitting sources with both criteria. It can be 
seen that the results are very close to each other. 

After the rejection of the sources with large deviations, 55 of the 7/zm fields (and 94 
of the 15/xm fields) were discarded using the following two criteria in combination with 
the preceding selection of sources: a field must have (I) at least 30 sources remaining and 
(2) values of Ejk s spanning more than 2 magnitudes, i.e. the maximum Cjk s should be 
larger than at least 3.2, since the minimum of Cjk s is set at 1.2 mag which agrees with 
the choice of RGB tip or early AGB stars (see succeeding paragraphs for explanation) . 
The first condition is required for the results to be statistically meaningful. The second 
condition means that the extinction at 7 and 15/im is large enough to be detectable. If we 
adopt our previous result that A7 ~ 0.2Ejk s , then Ejk s must be > 2, i.e. Ay > 12mag, 
and A7 > 0.4, just about twice the ISOGAL photometry rms accuracy at 7/im. 

These conditions retain about 73% of the 7fim fields and 57% of the 15^m fields. The 
resulting set of data covers 151 fields (including 35 FA and 116 FC fields) suitable for 
the 7/im determination and 126 fields (including 34 FB and 82 FC fields) for the 15/xm 
determination, based on the ISOGAL-2MASS PSC. For the ISOGAL-DENIS PSC, the 
retained data cover 139 fields (36 FA plus 103 FC fields) with 7/im photometry, and 116 
fields (42 FB plus 74 FC fields) with 15/xm photometry 

Two examples showing rejected sources are shown in Fig. [2J In this figure, the left 
panel shows the sources from the field FC-02992+00012, which represents a good case 
with small error in k^; the right panel shows the sources from the field FC+01694+00081, 
which represents a mediocre case with relatively large error in k.7. 

In principle, both the slope and the intercept, thus the color excess ratio, can be 
derived from the fitting. However, we choose to fix the intrinsic color index C^ s7 or 
C^ sl5 . There are two reasons for doing this. First, the intrinsic color indexes C^ s7 and 
Cks15 °f RGB tip stars or early AGB stars have sm all scatter, with a dispersion of about 
0.2mag comparable to the photometric uncertainty (|Jiang et al I2OO3I) . Second, for many 
fields, when Cjj s7 and Ck; s15 are not fixed, the fitting leads to unrealistic values of G^ s7 
and C^ s7 because of the lack of good quality low extinction sources. The adopted values of 
C^- s7 and C^- sl5 are 0.35 and 0.40 respectively. The choice is based on the calculation of a 
4000 K blackbody radiation, which gives a value in good agreement with the observations 



in th e ISOGAL field FC— 01863+00035 with particularly good quality data l|Jiang et al 
20031) . In addition, CV. is also fixed to be 1.2 mag, since the observations to the fie lds in 



Baade's window by 



Glass et al 



(1999) and the modelling bv lBertelli et a l. (1994) both 
showed that this index of such late-M type stars is 1.2±0.1. Therefore, the fitted line is 
forced to pass through the point [0.35,1.2] or [0.40,1.2] in the K s -[7] vs. J-K s or K s -[15] 
vs. J-Ks diagram, respectively. Only the slope of the fitted line, k7 or kis, is derived 
during the fitting. Even with the previously described restrictions, the fitting results still 
appeared poor in some cases, which may be attributed to some restriction not being 
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appropriate. For example, C^ s7 is fixed, which is not really true, since C^ s7 is likely to 
be slightly different from case to case. The last selection is based only on fitting error, 
without relation to any physical meaning. The final sample requires a relative error of 
ky less than 5 percent, i.e. a^/kj < 0.05, which takes into account the error induced in 
the fitting procedure. This restriction excluded 22 fields. For the 15/im case, the relative 
error of kis is constrained to be <10% in order to retain sufficient data, as 6 fields are 
excluded. Thus results with small fitting errors are derived for 129 fields and 120 fields, 
around 7 and 15/xm respectively. 

From kf, the ratio of the color excess E(Ks-[7]) to E(J-Ks) is derived. To obtain the 
extinction A7, one more step is needed as can be seen from equation (2) in Section 3.1. 
The value of A7 depends not only on the slope k7, but also on the near-infrared extinc- 
ti on values, At and A k s - The influence of near-infrared extinction values was discussed 



Jiang et al 



1 2003). Generally, the extinction law in the near-infr ared is accepted as 



universal with no variation from sightline to sightline l|Draine 



Nevertheless, the 



measurements o f the near-in f rared extinction s differ a l ittle f rom group to group. The 



values given by 



while the values given by 



( 2003) adopted the former while 




van der Hulstl l|1946|) are in good 



(RfcL) are slightly highe r. 



Indebetouw et al 



agreement 



Jiang et al 



( 2005) adopted the 



later. The averag e val ues of A7 expressed in Ak s are 0.34 when using the values by 



Glass et al 



1999) and 



van der Hulstl l|1946|) . while being 0.46 when using the values by 



Rieke fc Lebofskvl l| 19851) . For reason of consistency, we would prefer to use the Glass fig- 



ures. However, this group of values give rise to some negative A7 (7/zm) values, which is 
obviously unreasonable, and may mean the Glass values u nderestimate the ne ar-infrared 
extinctions. So instead of the previous choice, made by bv lJiang et alJ l|2003f) . the near- 
infrared values by R&L are used. Although the extinctions at 7/rni relative to Ay depend 
somewhat on the near-infrared extinctions, the value normalized to the extinction in the 
Ks band is less dependent on it. Therefore, we express the extinction at 7/xm in terms 
of the extinction in the Ks band. In Table ^ the results from different near-infrared 
extinction values are shown in detail. 



4. Results and Discussions 

4.1. Extinctions at 7\im 

For the 115 of all the selected fields observed by both DENIS and 2MASS, we have 
attested that the values derived for Aj/Ak s from the two sets of near-infrared data 
are consistent. They display the same average value with an mis difference equal to 
0.04. Only two or three fields displayed anomalous differences for A^/Ak s and these 
arise from a problem with the DENIS data. Thus, the two sets of data yield practically 
identical results. However, since 2MASS completely covers the ISOGAL-DENIS fields 
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Table 1. Average extinction values in the LW2 (7/im) and LW3 (15/zm) bands based on 
three different near-infrared extinctions 

Band A x /A Ka 

vdH 1 ISG 2 R & L 3 
J 2.816 2.876 2.509 

K s 1 1 1 

7 0.37 0.35 0.47 

15 4 0.30 0.28 0.40 

1 Values for J and Ks derived from van de Hulst (1946) curve 

2 Values for J and Ks derived from Glass (1999) (update of van de Hulst's values) 

3 Values for J and Ks derived by Rieke & Lebofsky (1989) for stars towards the Galactic center 

4 The average is derived only from the LW3 fields. See Section 4.3 for reference. 

and includes some northern ISOGAL fields not included by DENIS, the results derived 
from the ISOGAL-2MASS PSC represent the entire sample and only they are considered 
hereafter. 

From the above choices of the intrinsic color indexes C^ s7 and C^- sl5 as having the 
values 0.35 and 0.40 respectively, C7 15 is inferred to be 0.05. This corresponds to the 
case of stars with very little dust. However, the criterion C715 < 0.4 was set to exclude 
the effect of circumstellar dust on the brightness at 7/xm and 15/Ltm, on the color indexes 
Cks7 or Cks15, and thus on the result. The critical value of C715 affects the A7 result 
only slightly, while A15 changes a little more. When the critical value is changed from the 
present 0.4 to 0.6, the average of A7 is decreased by only 0.01, and the average of A15 is 
decreased by 0.08; when it is changed from 0.4 to 0.2, the average of A7 is increased by 
0.004, and the average of A 15 is increased by 0.09. Such tendencies are expected: radiation 
by circumstellar dust increases the 7/im and particularly the 15/im fluxes, which imitates 
a decrease in the corresponding extinctions. Thus, including more of the sources with 
circumstellar dust envelopes would underestimate the extinction. As the result on A7 is 
little affected, the value chosen appears to be small enough to exclude the effect of dusty 
envelopes. From the fact that A15 is affected when the critical value is changed, it can be 
seen that this criterion cannot completely remove the influence of dust shells. However, 
if this critical value is moved towards the blue, the number of sources and fields would 
decrease significantly and thus lower the statistical significance of the results. Therefore, 
the value of 0.4 is kept for both bands accordingly. 

The extinction values A7 are displayed in Tabled [the full table is given for the selected 
129 fields in the electronic version at CDS (http://vizier.u-strasbg.fr/viz-bin/VizieR. 
From left to right, the table includes the ISOGAL field name, filter ("2" stands for 
"LW2" and so on), the number of sources included in the final fitting, the slope k.7, the 
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Fig. 3. Distribution of all the A7 values as a whole (left panel) and in the different filters 
(right panel). 



1-sigma error of k 7 , and the extinction value at 7/j.m relative to Ak s - There are 106, 16 
and 7 fields with the filters LW2, LW6, and LW5, respectively. 

From Tableland Fig. [31 it can be seen that A7/AKS has a distribution which can 
be well approximated by a Gaussian function, with the central value of 0.47Ak s and 
a FWHM of about 0.07Ak s - The minimum and maximum are 0.33Ak s to 0.55Ak s re- 
spectively. Such variation includes the differences in filter wavelength ranges, statistical 
"noise" in determining k7, systematic errors and the possible variation of A7/AKS along 
various lines of sight. We will try to disentangle these various contributions in the two 
following sections. 



4.1.1. The central value of A7 



The average value of the distribution of A7 (c.f. Fi gEQ) we refer to as the central value is 



0.47Aks- This value is very consistent with that of 



Lutz et al 



1 1996) - an extinction of 



0.45Ak s at 6.8^m. It also agrees well with the value 0.48Ak s at both 5.7 urn and 7.8/xm , 



Indebetouw et al 



(2005). 



derived for the I — 42 deg sightlinc using the GLIMPSE data by 

As we mentioned in the Introduction, the ISOGAL observation used three different 
filters centered at about 7//m, viz. LW2, LW5 and LW6. To see if the extinction around 
7^im varies with wavelength, we have first determined the mean values of extinction in the 
three bands. They are 0.48, 0.42 and 0.43A Ks for LW2, LW5 and LW6, respectively. The 
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are shown are the average ones. The horizontal lines denote the spectral coverages of the 
ISOCAM and Spitzer IRAC filters. The vertical bars show the ranges of the values of 
A7 and A15. The dotted curve displays the silicate features calculated from Rosenthal's 
analytic formula. 



standard deviations of A7 in these three cases are 0.016, 0.020 and 0.019, significantly 
smaller than the variation ranges of A7 derived for the corresponding filters. The results 
are plotted in Fig EH and compared with t he values ba sed on Spitzer/GLIMPSE data by 



Indebetouw et al 



( 2005) and ISO data by 



Lutz (1999). It is seen that our results are in 



good agreement with others. 



4.1.2. Variation of A 7 

In order to discuss whether the spread in the values of A7/AKS exceeds that due to 
the data noise, and whether it provides evidence for extinction variations along differing 
lines of sight, we first focus on the LW2 data which are by far the most complete and 
homogeneous, with the best quality. The filters LW6 and LW5 were used only in the 
most difficult and crowded fields to reduce the effects of saturation by strong sources 
which are particularly numerous in star forming regions. The distribution of the values 
obtained with LW2 for A7 / Ak s is reasonably well fitted by a Gaussian distribution with 
a central value of 0.48 and an r.m.s. of 0.03 (see FigEJl - 
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Table [21 gives the value of the r.m.s. a^ 7 of the slope IC7 from the linear fitting of 
the relative reddening J-Ks/Ks-[7] (Fig. |2J). The r.m.s. values include mainly the noise 
appropriate to the various magnitudes, but also other kinds of variations such as uncer- 
tainty in the intrinsic colors, source variability, extinction variation along a given line 
of sight, etc. For all these 7/mi fields, the average of <Jki is 0.011, while IC7 ranges from 
0.30 to 0.44, with an unweighted average of 0.355 (weighted average 0.352) and standard 
deviation 0.025. Such a value (0.025) appears significantly, albeit somewhat marginally, 
larger than the r.m.s. ^0.011 of k-?. There is thus some evidence that A7/AKS may vary 
significantly, depending on the line of sight. 

The range [0.37Ak s , 0.55Ak s ] (o r even 0.33Ak s if one inc ludes LW6 fields) of A7 from 



our analysis covers that obtained by 



Indebetouw et al 



( 2005) from 0.37 at the I = 284 deg 



sightli ne to 0.48 at the / = 42 deg sightline at 7.8/im. The value of 0.40A Ks around 7^m 
from Lutz (jl999^) is also inside this range. 

It has been known for long time that the selective extinction in the visual and ultravio- 
let range varies from one sightline to another, depending on the interstellar environment 
through which the starlight passes. The representative parameter H=Av/(Ab — Ay) 
changes from values as low as 2.1 (toward very diffuse interstellar clouds) to values as 
large as 5.6-5.8 (toward very dense interstellar clouds). Such variation has been explained 
by different size distributions among the grains, i.e. the larger Ry distribution has sig- 
nificantly fewer small grains than the smaller Ry distribution, which is expected since 
generally there is relatively less extinction at short wavelengths (provided by small grains) 
for larger values of Ry . 

However, in the near-infrared, the available observa tions are consistent w ith a uni- 
form extinction law at wavelengths shortward of 3.5/im. 



Cardelli et al 



( 1989) concluded 



that the shape of extinction law for 0.7/im < A < 3.5^m was independent of Ry within 
the errors of measurements. It is further inferred from the constancy of the near-infrared 
extinction law that the size distribution of the largest particles are almost the same in all 

j — n 

directions (Mathis 1990). Following this opinion, one would initially expect no variation of 
the extinction law in mid-infrared, since the mid-infrared extinction must be attributed to 
larger grains. However, with the gra in models which have a size distribution of log- norma l 



form (We ingartner fc Draine 



2001) rather than classical power law form, 



Drain! (2003) 



(Figure 4) showed that the extinction law in mid-infrared could be much more dependent 
on Ry even though the law in near- infrared is rather independent of Ry. Similar to the 
situation at UV and visual wavelengths, the larger Ry value results in larger extinction 
in the 2.5-8/im range. At the wavelength 7/^m of our interest, the extinction changes 
by a factor of about two when Ry varies from 3.1 to 5.5. Thus in principle, the varia- 
tion of the selective extinction around 7/mi is understandable, considering the variation 



in the dust grain size distribution. On the observational side, 



Indebetouw et al 



(2005) 



found slight but significant differences in the wavelength dependence of the mid-infrared 
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extinction from two regions. One possible explanation for this discrepancy between the 
near-infrared constancy and the mid-infrared variation of the extinction law is that the 
extinction at 7/im is already affected by the silicate 9.7/im feature. While the continuum 
extinction (attributed to the graphite grains, which are the only source of extinction in 
the near-infrared wavelength range) does not vary from sightline to sightline, the silicate 
feature extinction may change in different environments. It is well accepted that there is 
variation in Ay/Arg.7. St udies of the 9.7/xm feat ure yielded Av/ATg.7 = 9 towards the 
Gal actic center direction I Roche fc AitkerJll98ril) and about 18.5 toward WC stars (see 



e.g. 



van der Hucht et al 



tio varies substantially (Whittet 



199 



3a 



In addition, within the Taurus molecular cloud, this ra- 



1988). Such variations in strength of the 9.7/im feature 
and its profile would surely affect the selective extinction in mid-infrared. Nevertheless, 
the wavelength limit of our 7/im filters remains less than 8.5/im and is thus only at the 
edge of the silicate 9.7/xm feat ure, where the relative extinction caused by silicate to the 



peak (A^/Ap ea k) is below O.l ljDraine 



1989). The influence of the silicate dust on the 



7/xm extinction should therefore not be very significant. This may be the reason why the 
variation of the 7/im extinction is not very large. 

In order to look for further evidence that the variation of the 7/Ltm extinction law is 
caused by variable silicate extinction, one could try to compare the extinctions deduced 
for the different filters used in the ISOGAL observation. Because the silicate feature 
centers at 9.7/im and starts from about 7.5/im, a small part of the LW2 filter and a 
major part of the LW6 filter's bandwidths overlap the silicate feature, while the narrow- 
band filter LW5 is completely outside it. In the right panel of Fig. [3J the distributions of 
A7 in three filters are shown. 

However, it is not certain that such a comparison is relevant because the number of 
LW6 and LW5 observations is relatively small, with data of poorer quality, and they cover 
atypical regions of the sky such as the Galactic Center neighborhood and star forming 
regions. There is no obvious difference between the LW5 and LW6 values; however, the 
number of LW5 values is too small to infer any real conclusion. The values of A7 from 
the LW2 and LW6 data both have approximately Gaussian distributions with similar 
widths. The most distant point in the A7 distribution of the LW6 observations is from 
the field FC+01694+00081. From Fig. Hit can be seen that the fitting for the extinction 
is of relatively poor quality. Indeed this field is located in the direction of the extremely 



young open cluster NGC 6611 ( B elikov et al 



1999), which may be the reason for the 



large scatter in the J-Ks/Ks-[7] plot and thus the poor quality of the fit. The difference 
between the LW2 and LW6 groups is apparent in the central values of the distributions. 
The LW6 group values are centered at 0.43Ak s and the LW2 group values center at 
0.48Aks, i-e. the extinction in the LW2 band is slightly larger than that in the LW6 
band. However, it seems very difficult when explaining the origin of such a difference, if 
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Fig. 5. Distribution of A7 with Galactic longitude (left panel) and latitude (right panel). 
Different symbols denote different filters: dot for LW2, triangle for LW5 and circle for 
LW6. 



real, to disentangle a possible wavelength effect from the fact that the LW6 observations 
are towards very special regions. 

4.2. Correlation with locations in the Galactic plane 

Shown in Fig.[5]are the distributions of A7 with Galactic latitude and longitude. Because 
the ISOGAL project targeted fields in the Galactic plane, the latitude range is small, 
and the distance from the midplane is always less than 1°. There is no systematic change 
of A7 with latitude seen in this latitude range. 

The longitude range of the fields we have retained covers about — 60°to +60°. There 
is no clear evidence that the extinction changes systematically with longitude within 
this large range. The measurements in the bulge direction are abundant and the results 
look rather scattered. The random variations of A7 may arise partly from "noisy" data. 
However, they are not incompatible with moderate local variations of the extinction 
caused by patchy distribution of the interstellar medium on small scales. 



4.3. Extinctions at 15fim 

Different from the 7/xm extinction, the 15/im extinction is probably dominated by the 
combined effect of spectral features of silicate dust with a minor contribution from con- 
tinuum extinction. If the continuum extinction followed the power law A\ oc A -1 ' 7 , it 
would amount to 0.04Ak s at 15/im. However, the silicate spectral feature s can bring about 



much larger extinction. For example, according to the analytic formula of Rosenth al et al 



1 2000), the extinctions at 15/xm caused by the 9.7/zm and 18/um features are <~ O.OAAk 
and ~ 0.30Ak respectively. If further the spectral response of the ISOCAM filter is taken 
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into account, the convolution of the spectral transmission profile of the LW3 and LW9 
filters with Rosenthal's formula yields 0.44Ak s and 0.39Ak s respectively. Thus the 15/xm 
extinction is mainly contributed by the silicate spectral feature extinction, with the 18/im 
feature being the major contributor. In addition, because of the variation of the silicate 
features in different environments, the Ibfim. extinction law is expected to vary as well. 

Because the flux at 15/zm is more affected by circumstellar dust than at 7/xm, exclud- 
ing the sources with relatively thick circumstellar envelopes becomes more important. 
For this reason, only the FC fields observed at both 7 and 15/Ltm with [7]-[15]<0.4 are 
chosen. Using this constraint, most of the sources with thick circumstellar envelopes are 
expected to be rejected. We adopted an intrinsic color index C^ sl5 of 0.40, because a 
4000 K (typical for a red giant star) blackbody radiation has this value and it is also 
the result d e duced from the previously studied ISOGAL field with good quality data by 
Jiang et al. I^f 



The final results on the 15/xm extinction derived from the FC fields are shown in 
Table [the full table is given for the selected 76 fields in the electronic version at 
CDS (http://vizier.u-strasbg.fr/viz-bin/VizieR and displayed as a histogram in Figure 
EJ Among the 76 FC fields for measuring A15 that remain, 69 were observed with the 
LW3 filter and 7 with the LW9 filter. For the 69 LW3 fields (i.e. observed with the 
LW3 filter), the average value of A 15 is 0.40 Ak s with a standard deviation of 0.07. 
This is remarkably consistent with the expected values from the superposition of the 
continuum extinction and the silicate spectral extinction as analyzed in th e beginning 
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of this section. This result is also in agreement with the value 0.41Ak s from 
(2003) for one line of sight. As cta 15 =0.037, and A15 ranges from 0.19Ak s to 0.54Ak s , 
which is larger than 3<r, it is possible to attribute part of the range to the field-to-field 
differences. This is consistent with the expectation that the silicate feature varies with 
the interstellar environment. However, the distribution of A 15 is asymmetrical (see Figure 
EJ); most of the results are within 3er. Taking into account the much larger rms of A 15 , 
the reality of the variation of A15 must be proved more rigorously. For the LW9 fields, 
the small number of fields makes the statistics less meaningful. Nevertheless, the average 
value of A15 for them is 0.14Ak s with a standard deviation of 0.09. Shown in Fig. 
are the distributions of A15 for the LW3, LW9 fields separately and together. As can be 
seen in Fig. the value for the LW9 fields is clearly smaller than for the LW3 fields 
and significantly smaller than the theoretical expectation. The difference in A15 between 
LW9 and LW3, as well as the difference between LW9 and the value from Rosenthal's 
analytic formula, is larger than the uncertainty of the fitting error. It may be concluded 
from the former that the extinction decreases from 14.3/im (the reference wavelength of 
LW3) to 14.9/im (the reference wavelength of LW9) and, from the latter, that Rosenthal's 
expression overestimates the extinction around this waveband. But we should be cautious 
about the extinction result at LW9, because: (1) The LW9 fields are relatively peculiar 
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Fig. 6. The histogram of the selective extinction at 15/im. The solid line shows the 
results from the fields observed in the LW3 filter, the long-dash line the LW9 filter, and 
the dotted line the sum of the two cases. 



regions with generally strong star forming activity so that the extinction may differ from 
diffuse interstellar regions. As Rosenthal's analysis is towards an area with an Ak s of only 
one magnitude, the difference may be attributable to the different environment. Further 
observations of some diffuse regions are needed. (2) The radiation from the circumstellar 
dust is strong around 15/im, and its effect can not be purged from our selection of the 
color index, as discussed in Section 4.1. 



5. Summary 

The ISOGAL project surveyed about 20 square degrees in the Galactic plane and detected 
numerous evolved stars with similar intrinsic infrared colors. Based on the point sources in 
the ISOGAL PSC catalog, we have built a cross-identified catalog between the ISOGAL 
and the 2MASS PSCs similar to the ISOGAL-DENIS catalog of Schuller et al. (2003). 
It includes about 60,000 sources. After selection from the mid-infrared and near-infrared 
colors and photometry quality flags, the extinctions around 7 and 15/im were derived 
towards about 129 fields. The central value of the 7/im extinction is about 0.47Ak s 



Rieke fc Lebofskv 



with a width of about 0.07Ak s (using the near-infrared extinctions of 
1 1985)). The average extinction shows some variation in the small range of wavelength 
defined by the LW2, LW5 and LW6 filters, with the value from the LW2 filter slightly 
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larger than from the LW5 and LW6 filters. In addition, there is a possibility that A7 varies 
along different lines of sight, falling within the range of [0.33, 0.55]Ak s , with the mean 
fitting error for any particular field being a bout 0.02. The valu es from H recombination 



lines toward the Galactic Center found by 



Lutz et al 



1 1996) fall within the range of 



our results, as does the recent result from the Spitzer GLIMPSE data at I = 42deg 



( Ind ebetouw et al 



2005) 



The extinction at 15/im was found to be about 0.4Ak s based on the data from fields 
observed both at 7 and 15/Ltm. Although this result suffers relatively large uncertainty 
due to scatter in the intrinsic color C^ sl5 , it is in reasonable agreement with theoretical 
expectation. 
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Table 2: Extinctions at 7/im from the selected ISOGAL fields 



Field Name 


Filter 


Number 


k 7 


0"fe 7 


A 7 


FA+00580+00076 


2 


155 


0.37 


0.012 


0.44 


FA+007f2+00068 


2 


80 


0.37 


0.010 


0.44 


FA+00737-00070 


2 


287 


0.34 


0.009 


0.48 



Table 3: Extinctions at 15/im from the selected ISOGAL fields 



Field Name 


Filter 


Number 


kis 


<7fc 1B 


Al 5 


FC-01904+00013 


3 


58 


0.41 


0.022 


0.38 


FC+00440-00009 


3 


36 


0.41 


0.028 


0.38 


FC+01514+00075 


3 


51 


0.40 


0.037 


0.40 



